Degree of randomness: numerical experiments for astrophysical 
signals 
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Abstract. - Astrophysical and cosmological signals such as the cosmic microwave background 
radiation, as observed, typically contain contributions of different components, and their statistical 
properties can be used to distinguish one from the other. A method developed originally by 
Kolmogorov is involved for the study of astrophysical signals of randomness of various degrees. 
Numerical performed experiments based on the universality of Kolmogorov distribution and using 
a single scaling of the ratio of stochastic to regular components, reveal basic features in the behavior 
of generated signals also in terms of a critical value for that ratio, thus enable the application of 
this technique for various observational datasets. 



Introduction. — Kolmogorov stochasticity parame- 
ter approach has been used to quantify the degree of ran- 
domness of sequences of number theory or dynamical sys- 
tems [iHS] . The stochasticity degree of the fractional parts 
of the arithmetical progressions has been analysed, also in 
comparison with geometrical progressions [5HZ] • 

The crucial aspect of the application of the method is 
the behavior of the stochasticity parameter. Arnold stud- 
ied the randomness for arithmetical progressions of the 
residues for the division by a real number and other se- 
quences, using the uniform or the Legendre-Chebyshev 
distributions [51|S] , mentioning also the problems for which 
no solutions are known yet. 

Along with this, physical applications of this descrip- 
tor of the degree of randomness have been also under- 
taken. Namely, astrophysical signals which contain several 
subsignals, i.e. contributions of different physical mecha- 
nisms, are of particular interest. Cosmic microwave back- 
ground (CMB) radiation, as observed, does contain be- 
sides the cosmological signal, also contributions of Galac- 
tic and other origins [Hin]. It appears that the different 
degree of randomness quantified by Kolmogorov distri- 
bution enables to distinguish the cosmological and non- 
cosmological signals (TOlIll] . The fractions of the random 
and correlated Gaussian components in the CMB over- 
all signal have been obtained using that method. It was 
also efficient for the identification in the CMB maps ob- 
tained by the Wilkinson Microwave Anisotropy Probe [12] 
of point sources [13j . i.e. of radio, as well as gamma-ray 



sources observed by the Fermi satellite. 

The properties of the CMB signal follow from primor- 
dial fluctuations which appear to be close to Gaussian 
ones. This remarkable empirical fact is fortunate for our 
approach, since it unambiguously defines the theoretical 
cumulative distribution function required to compute the 
stochasticity parameter. Although certain non-Gaussian 
features have been observed in the CMB datasets (see 
[T1HT7] ) and they are studied by various descriptors, the 
Gaussianity remains as one of the robust characteristics of 
CMB. However in other astrophysical signals, for example, 
for Lyman-alpha forest, the mentioned distribution func- 
tion cannot be so well defined and therefore the problem 
of statistics becomes more complex. In such conditions, 
the numerical experiments are the means to explore the 
behavior of the stochasticity parameter vs the generated 
random and regular components of the signal. Thus the 
universality of Kolmogorov's theorem is elaborated in the 
numerical experiments. 

Kolmogorov statistic. — First, let us briefly review 
the method. Consider {Xi,X2, ■ ■ ■ , Xn} independent val- 
ues of the same real-valued random variable X in growing 
order J^i < < • • • < X„ and let [Tl[3] 



F{x) = P{X < x} 



(1) 



be a cumulative distribution function of X . Their empir- 
ical distribution function Fn{x) is defined as 
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Fig. 1: The distribution of Kolmogorov's function <J> for the fuU 
sky cosmic microwave background temperature data obtained 
by Wilkinson Microwave Anisotropy Probe. 



Fn{X) = 

Kolmogorov's stochasticity parameter is 



0, X <xi 
k/n, Xk < X < Xk+i 

1, Xn<X 



A„ = Vn sup \Fn{x) - F{x) 



(2) 



(3) 



Kolmogorov theorem [T] states that for any continuous 



lim P{A„ < A} = $(A) 



where (E>(0) = 0, 



+ 00 



A > 



(4) 



(5) 



the convergence is uniform and Kolmogorov's distribution 
$ is independent on F. 

Since the stochasticity parameter itself is a random 
quantity its probable values are defined by the distribution 
function $, namely, the interval 0.3 < A„ < 2.4, where A 
provides information on the degree of randomness. The 
sky map of the degree of randomness for CMB is shown in 
Fig. 1: the map contains both the cosmological, as well as 
non-cosmological signals, e.g. the clearly visible Galactic 
disk (for details see [H]). Below we will study more tiny 
appearances of the Kolmogorov function vs the properties 
of the signals. 

Stochastic vs regular components. Using the 
Kolmogorov method we will study the properties of se- 
quences of type 



Zn = aXn + (1 - Ct)yn, 

where Xn arc random sequences and 

an (mod b) 
Vn = 7 ; 



(6) 



(7) 



are regular sequences, a and h are mutually fixed prime 
numbers, both sequences within the interval (0, 1) and 



have uniform distribution, a indicating the fraction of ran- 
dom and regular sequences. This representation is a more 
general case of random and regular sequences considered 
in [TTIIT5] . Other techniques for dealing with astrophysical 
systems with random (chaotic) and regular behavior are 
discussed in [T^ . 

When doing statistic with large number of sequences, 
each new sequence y„ is taken as the continuation of the 
previous one from the same arithmetical progression. 

Thus we have Zn with a distribution function 



F{X) = I 



2a(l-Q) • 
2aX-a^ 
2q;(1-q) ■ 

(1-XY_ 
2a{l-a) • 



0, X < 
Q<X <a 
a < X <l - a 

l-a<X <l 

1, X>1. 



(8) 



We will analyze the stochastic properties of z„ vs the 
parameter a varying between and 1 for different values 
of a and 6, i.e. corresponding to from purely stochastic to 
purely regular sequences. 

By fixing the values a and h, namely, a = 541, h = 1151, 
we generated 100 different sequences z„ for each value of 
a within and 1 by step 0.01, each sequence contain- 
ing 10000 elements. Then, each sequence is divided into 
50 subsequences and for each subsequence the parameter 
*&(A„)m is calculated (m runs through values 1, 50) and 
the empirical distribution function G'(<f>)m of these num- 
bers was constructed. In the case when the original se- 
quences are random, this distribution should be uniform 
according to Kolmogorov's theorem. Therefore, in gen- 
eral case we calculated of the functions G(<f>)m and 
Go(<i>) = $ to have an indicator for randomness. Thus 
one parameter is calculated per each of the 100 x 101 
sequences. Grouping 100 values per one value of a, we 
constructed mean and error values for ■ So, for each 
pair of a and b we get one plot: dependence of on a. 
Two examples of such dependence are given in FigO 

Figl2l shows the variation of vs a varying from to 
1, which corresponds to the gradual change of sequences 
Zn from regular to random. 

Then we study the behavior of this relation from param- 
eter a for different values of b. First we define a parameter 
A which equals to the difference of two values in the above 
plots: maximal value of x^ and minimal value in the range 
a € (0, amax), where amax is the position of the maximal 
value; e.g. it is obvious that for the left plot A = and 
for the right one A « 0.02. 

The next step is to fix b and for each value of primary 
a = 2, 6 calculate A. This is done for several values of 
b. 

FigElshows strict mirror symmetry in the dependence of 
A vs a, although no periodicity has been found by Fourier 
analysis. 

Then we proceeded via two types of analysis. In the 
first, we studied the values and the distribution of A and. 
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Fig. 2: The for the Kolmo gorov's function for the sequence 
Zn and of its subsequences vs a indicating the contribution of 
the random and regular parts for the pair of indicated input 
parameters. 
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Fig. 3: The dependence of A vs the parameter a. 
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Fig. 4: Sorted amplitudes of As from FigO 

in the second, the spacing between non-zero A and their 
distribution. Since obviously, most values of A are the null 
ones, we skipped them; also due to the mirror symmetry 
each A comes with its pair and therefore only one of each 
pair is taken into account. 

FiglHexhibits As sorted in growing order, i.e. on the ab- 
scissa axis we have the number of non-zero As from FiglSl 
The number of non-zero As appear to be proportional to 
b. 

Using the same values of As we calculated the frequency 
counts (FiglU, which showed the decrease of the number 
of non-zero As with the decrease of their values. 

The analysis of the spacings was done in the following 
way. We collected the distances of all neighbour non-zero 
As and constructed their frequency counts as shown in 
Figini smaller gaps are more common than larger ones. 
We also see the maximal gaps between non-zero As to be 
30 — 50, which do not change monotonically depending on 
a. 

Now let us follow the features of signals formed as sum 
of many fluctuations each having the same standard devi- 
ation. We generated two type of sequences, again of 10000 
elements each. 

The first one is given as 

1 ^ 

ai = -j= ^ Compact{xi, -1, 1), (9) 

k=l 

where = i/k, i.e. x'^ is compactified arithmetical se- 
quence within the interval (—1, 1), with step 1/fc. We call 
this a regular sequence. 
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Fig. 7: frequency counts for the set of random and regular 
sequences as compared with a Gaussian function. 



Fig. 5: The distribution of the amplitudes of A from Fig|3] 





Fig. 6: The frequency counts for the x-spacing of As. 



The second sequence is 
1 ^ 

hi = —= > Random(—l, 1), 

which we denote as a random one. 

Here we used the following notation: Compact{x,p^ q) = 
q + xmod{p—q) indicates multiples of (p—q) from x having 
the value within the range {p, q),p < q. The comparison of 
the sequences is done by means of Kolmogorov parameter 
with varying also the values of N. 

The results for 10000 random (generated by random 
number generator) and regular sequences each in FiglT] 
show the typical differences for when the number of the 
fluctuations vary from N = 1000 to 100000. The equal 
scales of shows that in both cases i.e. random and 
regular sequences, we deal with a Guassian distribution, in 
accordance with the Central Limit theorem. That theorem 
states that for large enough values of N, both sequences 
Oi and hi tend to Gaussian sequences with the same a 
and /Lt not depending on N. Along with that, however, 
we see differences in the Gaussians, namely, the standard 
deviations are larger for the regular case. 

Involving then the Kolmogorov function $, we see that 
although we have Gaussians both in random and regu- 
lar cases, the behavior of $ is quite different, as shown 
in FiglH i.e. it is close to a homogeneous function for 
random sequences and $ = 1 for regular ones. So, the 
Kolmogorov's function enables to distinguish the super- 
position of random and regular sequences, although both 
are tending to Gaussians. 

This analysis does not show any significant dependence 
on the length of the sequences, also dependence on the 
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Fig. 8: Kolmogorov function "I> for the sequences in Fig[7l 



number of the fluctuations within 1000 - 
weak as seen in FigE] 



100000 is rather 



Conclusions. — The performed numerical analysis 
aimed to reveal the behavior of the Kolmogorov distri- 
bution vs the randomness of generated signals. To de- 
scribe astrophysical datasets which contain both regu- 
lar and stochastic components, we considered sequences 
scaled by a single parameter a, the ratio of those compo- 
nents. 

Both qualitative signatures and their quantitative scal- 
ings have been observed at the numerical experiments. 
Namely, the appearance of a critical value of a has been 
shown, which defines the qualitative change in the behav- 
ior of the frequency count of the Kolmogorov distribu- 
tion: the monotonic decay is transformed to a function 
with extremum. Then, the dependence of these features 
vs the a reveals mirror properties in the amplitude and 
distribution of the frequency counts of the function $. 
Then, the behavior of the randomness for large number of 
subsignals has been also revealed, where the Kolmogorov 
function acts as an informative descriptor. Particularly, 
the descriptor is informative at large N limit when al- 
though the sum both of random and regular signals tends 
to a Gaussian in accordance to the Central Limit theorem, 
the set of regular subsignals is distinguished by a notably 
larger variance. In the case of the cosmic microwave back- 
gound, for example, the contributions to the overall de- 
tected signal, besides the cosmological one, typically are 
due to the Galactic synchrotron and dust emissions, star 
forming galaxies, point sources (quasars, blazars), instru- 
mental noise, etc. 

The behaviors obtained at the numerical experiments 



due to the universality of the approach will serve as in- 
dicators in the study of various astrophysical signals, in- 
cluding in revealing the contributions of correlated and 
random components. 

We are thankful to A.A.Kocharyan for many useful dis- 
cussions. 
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